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Angle-resolved photoemission spectroscopy (ARPES) is considered as the only experimental tool 
from which the momentum distribution of both the superconducting and pseudo-gap can be quan- 
titatively derived. The binding energy of the leading edge of the photoemission spectrum, usually 
called the leading edge gap (LEG), is the model-independent quantity which can be measured in 
the modern ARPES experiments with the very high accuracy — better than 1 meV. This, however, 
may be useless as long as the relation between the LEG and the real gap is unknown. We present 
a systematic study of the LEG as a function of a number of physical and experimental parameters. 
The absolute gap values which have been derived from the numerical simulation prove, for example 
that the nodal direction in the underdoped Bi-2212 in superconducting state is really the node — the 
gap is zero. The other consequences of the simulations are discussed. 

PACS numbers; 74.25.Jb, 74.72.Hs, 79.60.-i, 71.18.-l-y 



I. INTRODUCTION 

The question of symmetry of the superconducting gap 
is crucial to UBiJerstand the nature of superconductiv- 
ity in cuprates.tl Although not being phase sensitive, the 
angle-resolved photoemission spectroscopy (ARPES) has 
been considered as the only experimental tool from which 
both the superconducting and p|ttsCjudo-gap anisotropy 
could be quantitatively derived. EIltquIj The binding en- 
ergy of the leading edge of the photoemission snectrum, 
usually called the "leading edge gap" (LEG)Jj'EI is the 
model-independent quantity which can be measured in 
the modern AREJilS experiments with the highest, up to 
1 meV accuracy,!] and is, actually, the best quantity from 
which we can judge about the real gap values — in absence 
of a widely accepted model for the gap formation in these 
compounds, there is no direct way to extract the real gap 
value from the ARPES spectra. 

The LEG is determined as a lowest binding energy 
at which the energy distribution curve (EDC) reaches 
half of its maximum (here and thereafter we keep the 
same notation "LEG" for either gapped or non-gapped 
spectra). Therefore it is understandable, and usually ad- 
mitted, that the LEG should depend on any parameters 
that determine the EDC line shape and, unless the re- 
lation between the LEG and the real gap is known, can 
be considered only as a qualitative representation of the 
real gap. Moreover, the mentioned parameters (and EDC 
line shape) are substantially changing during any exper- 
iment, e.g. with temperature, along the Fermi surface 
(FS), etc. So, analysing the experimental data it is very 
important to distinguish between the artificial variations 
of the LEG and changes caused by the real gap in the 
electronic density of states. Despite the big importance 
of this question, there is no systematic study of the LEG 
available today. Hence, even the qualitative relation be- 
tween the LEG and the real gap can be doubtful. 

In this paper, by means of a numerical simulation we 
present a systematic investigation of the LEG as a func- 



tion of a number of physical and experimental param- 
eters which govern the ARPES spectra of the Bi-2212 
cuprates. 

Starting with the non-gapped case, we demonstrate 
that the dependencies of the LEG on the temperature 
and slight deviations in the momentum space (even 
within the momentum resolution window) from the FS 
(kp) result in strong variations of the LEG values which 
are comparable with the values of the superconducting 
and pseudo-gaps. We also show that at the low tem- 
peratures, choosing an appropriate criterion for the kp 
determination — the "maximum intensity" or "minimum 
gap locus" — the absolute value of the LEG is not very 
sensitive to the "physical" parameters of the model spec- 
tral function and therefore can be easily obtained. As 
an example we consider the experimental data of the 
LEG for jan underdoped Bi-2212 in the superconduct- 
ing statCjEI and prove that the nodal direction there is 
really the node — the gap is zero within the experimen- 
tal accuracy. Finally, switching to the gapped case, we 
show that the absence of a "cusp" of tke LEG vs FS 
angle dependencies (so called "U-shape"),Ll cannot be ex- 
plained neither by the interplay between the temperature 
and Fermi-function nor by other "artificial" parameters 
which may vary over the FS: momentum resolution, self- 
energy and band-structure. As an intermediate result, we 
present the tight-binding parameters for the considered 
compounds. 



II. MODEL DESCRIPTION 

We start with a model spectral function 



2 



which is an essential part of the photocurrent measured 
in the experiments 

/(k,cj,r, hv) oc [(Af(k, hiy)A(k,io,T)f{uj,T)] 

i?k«' R^ + B{uj,T). (2) 

Here M(k, hf) represents the dependence of the squared 
one-electron matrix element on in-plane electron mo- 
mentum k and excitation energy hv, f{Lu,T) = 
l/[exp(ci;/T) + 1] is the Fermi function. i?k and rep- 
resent the functions of the momentum and energy reso- 
lutions respectively. We include the energy resolution as 
a convolution with a Gaussian of the Ri^ full width at 
half maximum (FWHM) and the momentum resolution 
through a simple summation over a part of the detector 
area, 



A(»Rk = ^A, 



(3) 



<5k 



where (5k, in general case, is a 2D area with the R^ x Rk± 
dimensions (momentum resolutions parallel and perpen- 
dicular to the entrance slit of an analyser respectively) 
whose orientation in k-space depends on k and the ex- 
perimental geometry but, in the case of the big PS of 
superconducting cuprates, one can consider Jk as an one- 
dimensional cut ±i?fc(k)/2 perpendicular to the FS. 

For the momentum independent background we as- 
sume an empirical relation: B{uj,T) oc (1 -)- bijj'^)f{uj — 
Ab,T + Tb) with b = 1 eV-2, Ab = 5 meV and Tf, = 
90 K. The parameters At, and Tf, depend on doping and 
Ah also depends on temperatureO but we neglect this 
here because, as it is shown below, the typical signal- 
to-background ratios are large and the influence of the 
background on the LEG value is rather weak. 

For the imaginary part of the self energy we also use an 
empirical formula! S"(w,r) = {aoUjy + (0oT)^ which 
fits best the experimental data (see Refs. |9|,p^. In the 
vicinity of the Fermi-level (± 50 meV) the real part of the 
self-energy can be well approximated by its linear term 
E'(w, T) w -X{T)lo (A > 0) and the spectral function (0) 
can be rewritten in a renormalized form 



A{uj,k,T) 



|E"(^,T)| 



l + A(Lj-e(k))2+I]"(cj,T)2^ 



(4) 



with the renormalized imaginary part of the self-energy 



ao/(l + A), /3 = /3o/(l + 



A), and the renormalized dispersion e = e/(l -|- A). The 
last one, in a small vicinity of the Fermi energy Ep on 
the path perpendicular to FS, can be written using the 
renormalized Fermi velocity, e = vpk, which in the same 
way relates to the bare one: vp — uf/{^ + A). 

At last, when we take into account the bilayer splitting 
we include it as a simple superposition of the photocur- 
rent from bonding ("&") and antibonding ("a") bands: 

/ oc [MaA{ea) + MbA{eb)]f ®R],®R^ + B. (5) 



III. GAPLESS CASE 

First, we examine the LEG behaviour assuming that 
there is no real gap at all. It is natural to expect that the 
leading edge position of an EDC should depend on ev- 
ery physical parameter which forms the spectral function 
and the background as well as on energy and momentum 
resolutions. Among others, the temperature seems to be 
the most crucial here (see Ref. p^). 



A. Temperature dependence 

Fig. |l] shows the temperature dependence of the model 
photoemission spectra (a-c) defined by Eqs. ^ and ^ and 
the LEG values [d-f) derived therefrom. Here we use the 
following parameters: a = 1, (3 = 2, vp = eVA, R^ 
= 20 meV, and Rk ^ 0.07 A'^ that corresponds to 0.2° 
angular resolution for the nodal point at 21 eV excitation 
energy. Panel (a) represents the photocurrent (||) at fc = 
kp (kp-FiDGs) at different temperatures from 10 to 300 K 
(shown in color from blue to red with 10 K step). Panel 
(6) shows two extreme spectra but normalized to their 
maximum in order to illustrate the motion of the leading 
edge with temperature. The LEG(T) dependencies are 
shown on panel (c?) and (e) for different momentum and 
energy resolutions. 

The signal-to-background ratio (which can be esti- 
mated as a ratio of peak-to-tail photocurrent values) for 
the EDCs presented on panel (a) is chosen to be sim- 
ilar to typical photoemission data, and panel (e) also 
shows that the influence of the background in this case 
is negligible — the LEG(T) for the spectra without the 
background are shown as dashed lines. 

It can be seen that, although the LEG strongly de- 
pends on temperature and moves about 15 meV going 
from 10 to 300 K, which should be taken into account 
when the temperature dependence of the real gap is stud- 
ied, this dependence, for reasonable energy and momen- 
tum resolutions (R^ = 20 meV, Rk — 0.07 A~^: shown 
as a bold red curve in Fig. |lj d-f ), is quite monotonic 
(i.e. has no maxima on both itself and its derivative) and 
cannot imitate a rapid gap opening at a certain temper- 
ature. 

The situation is different if EDCs are taken from a k- 
point which does not exactly coincide with kp. Panels 
(c) and (/) of Fig. |^ illustrate this case. Panel (c) shows 
EDCs for k = kp + 0.01 A'^ for the same temperature 
range (10-300 K) and panel (/) shows the LEG(T) de- 
pendencies for different k — kp (from -0.01 to 0.01 A~^). 
One can see that the unoccupied part of the Brillouin 
zone (k > kp) is the most dangerous in this sense. If 
one steps away from fc^^^ in the unoccupied direction in 
0.01 A"i (about 0.3° for the nodal point at 21 eV ex- 
citation energy), the LEG changes about 30 meV going 
from 10 to 300 K but, what is most important, a kink ap- 
pears at about 80 K on the LEG vs T dependence. This 
kink can be easily misinterpreted as a gap opening, so we 
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FIG. 1: Temperature dependence of the model photoemission spectra with no gap included (a-c: temperatures from 10 to 
300 K shown in color from blue to red with 10 K step) and the leading edge gap (LEG) values derived therefrom for different 
momentum (rf) and energy (e) resolutions and positions in k (/). The lower curves on panels a and c represent the temperature 
dependent background. The red bold curves on panels d-f show the LEG for typical experimantal parameters Ru, — 20 meV, 
Rk = 0.07 k~\ k = kp. 



may conclude that the uncertainty in fci? determination 
(see discussion in Rcf. |l0| and below) is the most crucial 
parameter for the correct LEG evaluation. 

Regarding Fig. |l| (a) we also would like to note that due 
to the temperature dependent self-energy the area under 
each EDC is not constant with temperature, as one could 
expect according to the corresponding sum rule (e.g., see 
ReL |l|). 



the former and by the angular resolution perpendicular 
to the slit for the latter. As a result of this, the related 
momentum resolution can change by a factor of two but, 
as one can see from Fig. (d), this will move the lead- 
ing edge about 1 meV only (even for the "worst" case of 
10 K). 



Band structure 



B. Resolution 

The dependence of the LEG on resolutions can be es- 
timated from Fig. |l| (rf) and (e). The dependence on 
energy resolution (Fig. ^ (e)) is stronger but it is not a 
problem because we believe that normally it remains con- 
stant during the experiment and the absolute value can 
be easily taken into account. The momentum resolution 
can change during the experiment. For example, for the 
experimental geometry described in Ref. 10, the actual 
momentum resolution changes going from the r-(7r, tt) 
to (tt, 7r)-(7r, 0) FS crossings: it is defined by the angu- 
lar resolution along the entrance slit of the analyser for 



In contrast to the temperature, the other parameters 
that enter the fci?-EDCs are expected to depend on the in- 
plane momentum of the electrons. One can divide these 
parameters into "structural" , which come from the bare 
electronic structure, and "physical" , which come from 
interactions and reveal themselves as a self-energy. 

The structural parameters of interest here are the bare 
Fermi velocity and the bilayer splitting. Fig. ^ represents 
the result of the tight-binding fit of the experimental data 
published in Refs. and |ll[ The precisely determined 
FS of Bi-2212 (at least its bonding sheet, see Ref. |ll|) at 
room temperature (300 K) can be well described within 
the four-band model introduced in Ref. Ol Within this 
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FIG. 2: The result of the tight-binding fit to the bonding (red solid lines) and antibonding (blue dashed lines) CuO-bilayer 
bands of the overdoped (OD 69K) Bi-2212: left panel shows the corresponding Fermi surfaces (the boundary of the first Brillouin 
zone marked by dotted square), right panel shows the 'bare' dispersion along the F-Y-M-F path (shown in the left panel as a 
triangle), and inset zooms in the (7r,0) region, from (7r,7r/6) to (tt,— 7r/6). 



model, the simplified bare dispersion 

ea,b = Ae — 2t{cos kx + COS ky) 

+ At' cos kx cos ky — 2t"(cos 2kx + cos 2ky) 

± t ±{cos kx — cos ky)'^ /A, (6) 

where second and third nearest-neighbor intra-plane hop- 
ping integrals {t' and t") are provided by the "Cu s" or- 
bital and interlayer hopping is described by the t± (the 
bilayer splitting, the energy distance between—bonding 
and antibonding CuO-bilayer bands, is 2ij_)E3'EJ 

To fit the data we use the following procedure. First, 
the relative parameters (relative to t) are determined 
from the room temperature FS maps of Ref. then the 
bilayer splitting is taken into account [the value tj_/t is 
determined from the (tt, tt) — (7r,0) FS crossing]. For ex- 
ample, for the overdoped samplc-iOD, Tc — 69 K), using 
the assumption that t" « t'/2,t3 these parameters are: 
t'/t = 0.23, t"/t = 0.11, tj_/t = 0.21, Ae/t = 1.08. From 
here one can estimate the variation of the bare Fermi 
velocity up around the FS as a ratio of its value in the 
antinodal 'A' point (FS crossing along (tt, tt) — (tt, 0) di- 
rection) to the nodal ((0,0) — (7r,7r) crossing) one. This 
ratio, up{A)/uf{^), is 0.71 for the bonding band and 
0.25 for the antibondig band (for OD 69 K). We stress 
that these values are determined by the FS shape, as pre- 
sented on the left panel of Fig. ||, and do not depend on 
the energy scale. 

In order to estimate the "scale" t one needs to know 
some absolute value of the bare band. We can esti- 
mate the 'bare' Fermi velocity as 4.0 eVA, subtracting 
the real part of the self-energy, which has been deter- 
mined using the Kramers-Kronig transform of its imag- 
inary part, from the renormalized experimental disper- 
sion. The tight-binding parameters which have been cal- 
culated for two samples with this value of wf, assuming 



TABLE I: Tight-binding parameters of the CuO conducting 



band of Bi- 


2212. 








Sample 


t (eV) 


t' (eV) 


t" (eV) tx (eV) 


Ae (eV) 


OD 69 K 
UD 77 K 


0.40 
0.39 


0.090 
0.078 


0.045 0.082 
0.039 0.082 


0.43 
0.29 



the constancy of the bilayer splittingJ13 are presented in 
Table ||. The right panel of Fig. |^ shows the obtained bare 
dispersion along the F-Y-M-F path (see the left panel). 
The inset in Fig. |^ zooms in the Y-M-X region (exactly 
the (7r,7r/6) — (tt, — 7r/6) path). This gives the bilayer 
sphtting in (7r,0)-point of 160 meV and the absolute bare 
Fermi velocities: uf(A, bonding) = 2.8 eVA, u_f(A, an- 
tibonding) = 1.0 eVA, while uf(N) = 4.0 eVA. 

One can expect that the renormalized Fermi velocity 
vp = uf/{1 + X) varies even more strongly than the 
bare one going from the node to th e . an tinode due to a 
change of the 'coupling strength' A.E^Ea So, going back 
to the simulation procedure, we examine an effect of 
variation on the leading edge position in the range from 
0.5 to 5 eVA. The results are represented in Fig. a {a,b). 

From Fig. ^ (a) one can see that the influence of wf 
on LEG increases with the decreasing temperature but 
remains rather weak - about 4 meV change at 10 K over 
the whole examined range (0.5-5 eVA) or in the more 
reasonable range from 1.5 to 3.5 A (shadowed on panel 
b) inferred from the experiment,0 this variation is only 2 
meV for k = kp. Fig. ^ (6) also shows that if one follows 
the minimum gap locus (MGL) criteriaEl — choosing the 
location in k for which the LEG(A:) along a FS cut has 
a minimum — this change is even less but, if k deviates 
more {\k — kp\ > Rk/2 = 0.035 A^^), the change rapidly 
increases. 

The influence of the bilayer splitting on the leading 
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FIG. 3: Dependence of the leading edge position on the band structure parameters: (a, b) renormalized Fermi velocity, 
vf = uf / [I + A); and (c,rf) bilayer splitting, Ae = Sa ~ St, = Ma/Mi, (see Eq. and k = feF(bonding). 



edge is shown in Fig. || {c,d). Here Ae = Sa — Sb is the 
value of the splitting, and fi — Ma/Mb (see Eq. |^) and 
k = kp of the bonding band. One can see that at 120 
K the influence of other FS sheet (antibonding in this 
case) can produce some minor variations (within about 2 
meV) in the LEG around the FS but at low temperature 
(30 K) these variations are negligible. 



D. Self-energy 

The influence of the self-energy parameters on the lead- 
ing edge position is presented in Fig. ^. Although, due 
to the presence of the bilayer splitting, there is no sys- 
tematic information currently available about how these 
parameters change over the FS, one can expect that in 
the antinodal point the coupling strength is larger than 
in the nodal one and, therefore, the self-energy parame- 
ters a and (3 (see above) can be a few timearbigger — two 
times seems to be a reasonable estimation.Bll3 

Fig. ^ (a) shows that the variation of a even in a much 
wider range (from 0.2 to 4) moves the leading edge within 
about 1 meV only, if k — kp. Small deviation from the 
kp noticeably increases the amplitude of the LEG(q!) de- 
pendencies (Fig. 1^ (6)). 

The dependencies of the LEG on /3 within the same 
range, being more temperature dependent, are similar in 
amplitude (Fig. | {c,d)). 



IV. GAPPED CASE 

A. Minimum gap locus 

Among the kp determination criteria, the "minimum 
gap locus" (MGL) intuitively seems to be the most suit- 
able in a gapped state. On the other hand, as it is shown 
above, the MGL, if applied to the kp determination, gives 
overestimated values. What is important here is that the 
/c-locations which come from the MGL-criterion, kMGL, 



are in the "dangerous" region, where a strong depen- 
dence of the LEG on temperature and other parameters 
is expected. We apply the simulation to check at which 
conditions the MGL method can be used. 

Fig. H shows the dependencies of the leading edge posi- 
tion on k along a cut perpendicular to the FS for temper- 
atures from 10 to 300 K (from top to bottom) with 10 K 
steps. We use parameters typical for the F-X crossing: a 
^l,f3 = 2,vp^2 eVA, = 20 meV, Rk = 0.07 A-^. 
Fig. |5| (a) represents the gapless case and Fig. || (6) shows 
the same but taking into account a gap which opens at 
Tc = 90 K. We used the BCS modified spectral function 



,2 V" 



A (X 



(7) 



where = ^4 A(T)2, ul^\{l + Sk/Ek), and the 
gap function for T < Tc is approximated as 



A(T) = Ao 



1.74W1 



(8) 



with Ao = 1.76 Tc (temperature in energy units) and 
n = -4. 

The gapped region is shadowed in Fig. || (6) with the 
blue color, and the gapped LEG(fc, 10 K) dependency is 
also shown in Fig. |^ (a) as a dashed curve to compare 
with the gapless case. 

From the results presented in Fig. || one can make 
the following conclusions. (1) The value of k at which 
the leading edge position reaches a minimum, kMGLj is 
strongly temperature dependent — thus it is important to 
realize that the MGL spectra measured at different tem- 
peratures correspond to rather different /c-points. (2) 
The difference between kMGL and kp is improperly big 
at high temperatures (above Tc) and remains quite big 
at lower temperatures for both gapped and non-gapped 
cases. Using a criterion \kMGL — kp\ < Rk/2 one can con- 
clude that the MGL method cannot be used for precise 
kp determination for T > 40 K, at least in the Bi-2212 
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case. (3) If one stays at a certain /c-point, measuring the 
LEG as a function of temperature, then the region where 
the real gap opening can be distinguished from the ar- 
tificial one can be estimated as fc < + i?fc/2. In this 
case the maximiim intensity methodJ13 which underesti- 
mates the kpf^ seems to be the most suitable. (4) At 
low temperatures (T < 40K) the MGL method is quite 
good as a method of the LEG determination — the exper- 
imental uncertainty in Umgl is big here because of very 
flat LEG(/c) dependence but, as another consequence of 
this, the LEG uncertainty should be small. 



B. LEG in the case of real gaps 

Now we consider how the real gap can influence the 
LEG value. Fig. ^ represents the results. Panels (a) 
and (d) correspond to kp-spectia and panels (6) and 
(e) show the LEG of the spectra from k = kMGL- 
The kMGL{^) dependencies at different temperatures are 



FIG. 5: The dependencies of the leading edge position on k 
along cut perpendicular to the FS for temperatures from 10 
to 300 K (from top to bottom) with the 10 K step for typical 
for Bi-2212 parameters {a = 1, {3 = 2, vf = 2 eVA, = 20 
meV, Rk ~ 0.07 A~^): panel (a) represents the gapless case, 
in panel (b) the BCS-hke gap opens at Tc = 90 K (the gapped 
region shadowed with the blue color). 



shown in panel (c). Note, that all values are shown as 
function of A(r) rather than Aq. 

If we accept the use of Eq. |^ for the gapped spectral 
function the question how the imaginary part of the self- 
energy depends on the gap is unclear — in other words, 
the I]"(A) function is considered to be "more model de- 
pendent" than the spectral function A{A). The LEG(A) 
dependencies presented in Fig. ||(a,6) have been obtained 
with a A-independent self-energy. In Fig. ^ {d,e) we 
model the influence of the gap on the imaginary part 
of the self energy as A,r) = y/{agJy + {pxy, 
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FIG. 6: The leading edge gap (LEG) and the minimum gap locus (MGL) vs real gap value A(r) for the models with the 
A-independent (a-c) and BCS-like {d,e) self-energies. 



where 

ifM>A; 

which gives the BCS asymptotic oc oj^ at uj ^ 

and zero teniperature. 

One can see that there is some "rounding" (deviation 
from hnear) of the LEG(A) curves at low A which is 
explained by cutting finite width EDCs by the Fermi- 
function. These "rounding" increases with temperature 
but it is possible to conclude that for low temperatures 
(T < 30 K) it remains rather small (less than 2 meV) for 
all fc-points and T,"{uj) models considered. Tiiis means, 
for example, that the U-shape of the LEG((/9)Ey observed 
in Ref. |^ cannot be explained by such a rounding of a 
sharp cusp expected in a simplest d-wave gap function. 

The comparison of Fig. ^ {a,b) and {d,e) respectively 
shows that at higher A the LEG(A) dependencies be- 
come model dependent. Then, from the "more physi- 
cal" model of T,"{uj), determined by Eq. ^, we can get a 
rough estimation of the coefhcient between the LEG and 
A: dLEG(A)/dA « 0.5. 

V. DISCUSSION 

With the above simulations we mostly focused on the 
gapless case — we examined the leading edge position of 
the photoemission spectra in case when there is no gap 
present in the electron density of states. 

It has appeared that the leading edge position of such 
a gapless spectrum depends on temperature, resolutions, 
band structure and self-energy parameters, which is not 
surprising of course. A surprising result is that under 
certain circumstances these dependencies can be quite 



complicated (exhibit non-monotonic rate) that can be 
misinterpreted as a physical transition. 

These "circumstances" are a wrong position in mo- 
mentum space, i.e. when fc-location of the spectrum does 
not coincide with kp but is uncertain. If the parameters 
mentioned above are well determined experimentally and 
can be taken into account with the described simulation, 
the uncertainty in momentum cannot be dealt with in 
the same way, by definition. More precisely, this uncer- 
tainty is the uncertainty in k — kp and, therefore, consists 
of the two: fc-uncertainty, which is determined by the k 
step of the measurements, and fci?-uncertainty which is 
related to a, kp determination procedure. The first one 
can be reduced by refining the k step of the measure- 
ments, which should be significantly less that the mo- 
mentum resolution. The second problem is well known 
and heavily discussed in the past but we would like to 
address this question a bit more. 

One way to define an uncertainty of the method is 
to model how far the directly determined value of kp 
deviates from the real one — for example, in the "max- 
imum intensity" method, how the value of k, at which 
the photocurrent intensity exhibits a maximum, devi- 
ates from the real kp. It has been shown that for high 
temperatures the maximum intensity method, especially 
in combination with the normalization to "highest bind- 
ing energy", is much more precise than others, sucji-aa 
"Vn(k)", "minimum gap locus" or "AT" methodsJ13'til 
but modelling this, one can find that it is difficult to sug- 
gest such a universal procedure of the kp determination 
which always gives the best result. All of the methods 
give some deviations and then one can reason that it is 
possible to determine the real kp taking these deviations 
into account turning the uncertainty to zero. Then we 
come to another uncertainty: how much these deviations 
depend on the model of the spectral function we use. 
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This is an interesting and complicated question which 
can be subject of a separate study. Here we apply the 
word "uncertainty" of the method in the first sense and, 
in this respect, the maximum intensity method appeares 
to be the most accurate. 

What is the real benefit one can extract from the re- 
sults obtained or what do these results contribute to the 
gap problem? There are two kinds of gap measurements 
which are considered to be "model independent" : the gap 
as a function of temperature, LEG(r), and momentum, 
LEG(k), which is usually presented as a function of FS 
angle, LEG((p). In both cases the obtained "gap" is rel- 
ative: in the first case it is relative to the LEG at the 
highest temperature, in the second case it is relative to 
the LEG in nodal direction. 

The presented simulation shows that even when a real 
gap is zero the LEG is strongy temperature dependent 
(possibly in a complex way), and the amplitude of this de- 
pendency is comparable with the estimated real gap val- 
ues in the superconducting cuprates. On the other hand, 
the momentum dependence of such a gapless LEG at low 
temperatures is shown to be rather weak for all possible 
variations of physical, structural and experimental pa- 
rameters. But what is most important is that knowing 
these parameters one can estimate the absolute leading 
edge position and, therefore, answer the question whether 
there is any gap in a given spectrum or not. 

For example, the LEG at 40 K and with the other pa- 
rameters used in Ref. R for the nodal direction is equal 
to —7 meV (see Fig. ^ d-f) which precisely coincides 
with the value obtained in the cited paper. This fact al- 
lowed us to make a firm conclusion that the nodal direc- 
tion in an underdopcd Bi-2212 in superconducting state 
is really the node — the gap is zero — which is in favor of 
d-symmetry of both superconducting and pseudo-gaps. 

Another conclusion comes from Fig. |: the "U- 
shape" of the LEG(i^) dependencies, wkich have been 
recenly observed for underdoped samples,Ll cannot be ex- 
plained neither by the interplay between the tempera- 
ture and Fermi-function nor by other "artificial" param- 
eters which may vary over the FS: momentum resolu- 
tion, Fermi velocity, self-energy parameters. Considering 
the absence of the gap in the node as an evidence for 
d-wave symmetry of the order parameter in Bi-2212, this 



U-shape can be a consequence of higher harmonics of the 
gap due to the long rarige of the pairing interaction, as 
it was suggested earlier ,0 or due to an interplay between 
the superconducting and pseudo-gaps. We want to note 
that this conclusion is valid for low temperatures only 
(superconducting state for Bi-2212), for higher temper- 
atures (e.g. for pseudo-gap state) the influence of the 
artificial parameters can be essential. 



VI. SUMMARY 

In conclusion, a systematic study of the LEG is pre- 
sented. Even in the absence of a real gap, the lead- 
ing edge position of the ARPES spectrum depends on 
many factors: energy and momentum resolutions, Fermi- 
velocity, the self-energy, but mostly on temperature. At 
some circumstances, when the spectra are measured away 
from kp (if k — kp > 0.005 A^^ for typical parameters of 
an experiment on Bi-2212), the temperature dependence 
of the LEG can be quite complicated and exhibit non- 
monotonic rate that can be misinterpreted as a physical 
transition. At low temperatures the dependence of the 
LEG on other parameters is rather weak and cannot be 
responsible for observed in experiments the LEG varia- 
tion over the FS. 

The absolute gap values which have been derived from 
the presented simulation prove that the nodal direction 
in the underdoped Bi-2212 in the superconducting state 
is really the node (the gap is zero) but also that the "U- 
shape" of the LEG vs Fermi surface angle, observed for 
underdoped samples, is not artificial. This implies the d- 
wave symmetry of the superconducting and pseudo-gaps 
but with essential contributions from the higher harmon- 
ics of one or both gaps which can be a result of the long 
range pairing interaction. 
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